Ground-state properties of two-dimensional dimerized Heisenberg models 
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The purpose of this paper is to investigate the ground-state properties of two-dimensional 
Heisenberg models on a square lattice with a given dimerization. Our aim is threefold: First, 
we want to investigate the dimensional transition from two to one dimension for three models 
consisting of weakly coupled chains for large dimerizations. Simple scaling arguments show that 
the interchain couphng is always relevant. The ground states of two of these models therefore 
have one-dimensional nature only at the decoupling point. The third considered model is more 
complicated, because it contains additional relevant intrachain couplings leading to a gap as 
shown by scaling arguments and numerical investigations. Second, we investigate at which 
point the dimerization destroys the Neel ordered ground state of the isotropic model. Within a 
mapping to a nonlinear sigma-model and linear spinwave theory (LSWT) we conclude that the 
stability of the Neel ordered state depends on the microscopic details of the model. Third, the 
considered models also can be regarded as effective models for a spin system with spin-phonon 
coupling. This leads to the question if a spin-Peierls transition, i.e. a gain of total energy due 
to lattice distortion, is possible. LSWT shows that such a transition is possible under certain 
conditions leading to a coexistence of long-range order and spin-Peierls dimerization. We also 
find that the gain of magnetic energy is largest for a stair-like distortion of the lattice. 
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I. INTRODUCTION 

Since many years there is considerable interest, both 
experimentally and theoretically, in the subject of 
low-dimensional quantum spin systems, because their 
properties are strongly affected by quantum fluctuations. 
The generic model to theoretically study such systems is 
the well-known Heisenberg model. 

In one spatial dimension the model with nearest 
neighbor exchange of spin-1/2 objects, known to be 
exactly solvable by Bethe ansatzjll shows an algebraic 
decay of its correlation functions at zero temperature 
and constitutes thereby a qviantum critical system. 
For arbritary spin s, Haldand3 has mapped the spin 
chain onto a nonlinear sigma-model (nlcr-model) with a 
topological term for half-integer s and without such a 
term for integer s. From this result he conjectured that 
half-integer spin chains are critical whereas integer spin 
chains have a gap, a scenario that is well established 
by now. Another interesting aspect of the spin-1/2 
Heisenberg chain is its instability towards-, a structural 
transition known as spin-Peierls transitionJa 
Much less is known for the isotropic two-dimensional 
Heisenberg antiferromagnet with nearest neighbor 
exchange on a square lattice. Contrary to one dimension 
no exact solution is available in any limit. LSWT, which 
does not work in one dimension because of infrared 
divergencies, is applicable and predicts a Neel ordered 
ground state for the spin-1/2 case, but with a magnetic, 
moment reduced to nearly 50% of its classical value.cl 
This result is also supported, quaJitatively and also 
quantitatively, hv numerical work.ErQ For s > 1 Dyson, 
Licb and SimonQ proved a theorem, which shows that 



the ground state is Neel ordered. Like in one dimension 
it is also possible to map the system onto a nlcr-model. 
From pa renormalization-group (RG) treatment it is 
knownB that this model exhibits in (2-|-l)-dimensions 
a nontrivial critical point gc, which separates a phase 
with Neel-like long-range order (g < gc) from a quantum 
disordered phase {g > gc) at T — 0. It has been shown 
that there is an excellent agreement between theoretical 
results for this model and experimental measurements 
on La2Cv,04 in the low-temperature regime ii g < gc is 
assumed.B An interesting problem has been the question, 
if there is a topological term also in|-t.wo dimensions, 
which was finally answered by HaldanejlJ who concluded 
that such a term is always absent if the order parameter 
field is smooth on the scale of the lattice spacing. 
However, there are tunneling events, which are crucial 
for an understanding of the disordered phase. til 

In this work we want to consider two-dimensional 
Heisenberg models on a square lattice with a given 
alternation of the coupling between nearest neighbor 
spins. In each spatial direction the coupling should 
be changing from bond to bond between J(l -I- i5) and 
J(l — 6) with J > and S e [0, 1] so that the couphng is 
always antiferromagnetic. There are three topologically 
different possibilities (cf. figures |^, ^ ^) for arranging 
such " dimerized chains" on a square lattice if periodicity 
in each spatial direction is assumed. These systems are 
described by the following Hamiltonian 



H ^ jJ2[l + i-iy(+^^5)s,,,S,+,,, (1) 

i 

+ J^(l + (-l)^(+')<5)s,,,5„+i, 
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FIG. 1: (STAIR-MODEL) A thick solid line indicates a strong 
bond J(l + 5) and a dashed line a weak bond with strength 
J(l — 5). This distortion of the lattice is caused by one 
transversal phonon with wave vector (7r,7r). 



FIG. 4: ( Ja;/ Jy-MODEL) Two longitudinal phonons, one with 
wave vector {—q, 0) the other with (0, +q), where q is infinites- 
imal {q = 2TT/VNa). 



FIG. 2: (MEANDER-MODEL) Three phonons, a longitudi- 
nal and a transversal (tt, 7r)-phonon together with a longitu- 
dinal (0, 7r)-phonon. 



FIG. 3: (PLAQUETTE-MODEL) Two longitudinal phonons, 
one with wave vector (tt, 0), the other with (0, tt). 



where Sij denotes the spin operator acting on the lattice 
site Choosing both exponents equal to i + j leads 

to the model shown in figure |l|, whereas setting the first 
one equal to i and the second one to i + j or vice versa 
leads to the model shown in figure |[ The third consid- 
ered model (see fig ^ is described by (|l|) with the first 
exponent set to i and the second set to j. 
The models in figure [l] and fi eure 2, which we want to 
refer to as stair-model and meander-model, decouple into 
spin chains at (5 = 1. This means that there is a transi- 
tion from two to one dimension depending on the value 
of the dimerization S. Because a model with a coupling 
Jx in x-direction and a coupling Jy in y-direction (see 
figure^) is the simplest model showing such a transition, 
we want to reexamine this model although it has been 
studied intensively before. The model in figure |^ is dif- 
ferent from the other three models, because it decouples 



into plaquettes consisting of four spins for ^ = 1. 
The dimerized models also can be regarded as effective 
models for a spin system with a spin-phonon coupling 
treated adiabatically. In such a system the exchange 
coupling J between nearest neighbors depends linearly 
on their distance. This is expressed through the defor- 
mation parameter of the horizontal bond extending from 
site (j,j) to the right 



) = J^>'k,s{{u^J{k,s)) - (Mi+ij(A;,s)))|| (2) 

k,s 



and the deformation parameter of the vertical bond 
S'"{i,j) being defined analogously. Here Xk s is the mi- 
croscopic spin-phonon coupling constant and Uij{k,s) 
the local displacement of the atom at the position 
with respect to the phonon wave vector k and branch 
s. In the considered models 5'' and S"" are equal to 
±6. Note that the lattice distortion is static due to the 
mean-field approach leading to a classical elastic energy. 

In section^ we study the models, which decouple into 
chains at (5 = 1. We address the question whether there 
is a transition from two- to one-dimensional behaviour 
already at a J smaller than 1 on the basis of some scaling 
arguments. Because the situation in the meander-model 
is much more complicated than in the other models due 
to third-nearest neighbor (NNNN) couplings, we show in 
section |l| numerical results from density-matrix renor- 
malization group (DMRG) and transfer-matrix DMRG 
(TMRG) to confirm the conclusions drawn from the sim- 
ple scaling arguments in section In section IV, we 
show that it is possible to map all models onto a 2-|-l- 
dimensional nlcr-model. We use the known RG results for 
this model to discuss the magnetic properties of the spin 
models. In section 0, LSWT is apphed to support the 
picture from the RG arguments and to determine a value 
for the critical dimerization Sc, where the magnetic order 
vanishes. Also from this, we get a condition for the spin- 
Peierls transition and an answer to the question, which 
structure is energetically prefered. The sections and 
^ consist of separate subsections for each model. How- 
ever, in we have interchanged the order of the first two 
subsections in comparison with |^ to keep calculations as 
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simple as possible. In section VI, we discuss our results 
and give some conclusions. 



II. SCALING ARGUMENTS 

At 6 = 1 — e with e ^ 1, the models in figure |l|, ||, 
^ consist of weakly coupled Heisenberg chains. Such a 
Heisenberg chain with s = 1/2 is a critical system and 
the additional weak inter- and intrachain couplings are 
small perturbations of this critical system. If we pick two 
such chains from each model with the corresponding in- 
terchain coupling proportional to e, we can determine the 
relevance of the perturbation by calculating the energy- 
energy correlation function and from this the scaling di- 
mension of the perturbation operator. Let us start with 
the simplest case, the Ja;/Jy-model, where Jy — J(l — S) 
and Jx — J(l + S). The Hamiltonian of the weak inter- 
chain coupling is given by 



there seems to be a gap for an even number of coupled 
chains, whereas a system with an odd number of chains 
is gapless. However, it is not possible to determine from 
scaling arguments if there is a gap or not for an infi- 
nite number of arbitrarily weakly coupled chains. Nev- 
ertheless, the relevance of the interchain coupling clearly 
shows that the system scales away from decoupled chains 
and therefore even at very large, however not perfect 
dimerization does not behave like decoupled chains. We 
conclude that the ground state of this model has twp^ 
dimensional nature if Jj, ^ as has been stated before.til 
Also we can pick two chains out of the stair-model and af- 
ter smoothing the chains, we get the configuration shown 
in figure ||. Here the interchain coupling is described by 



H ijy ^ ^ S J, S ^ , 



(3) 



FIG. 5: Two chains from the stair-model with the correspond- 
ing interchain coupling 



where the upper index labels the two chains. The energy- 
energy correlation function of this perturbation can be 
calculated as follows 



the Hamiltonian 



a, 13 



Jy 2_^W0 /Wo ^7- 



"^y / a Q \2 '^y ( 



(4) 



where Ur — JyS].Sl and a, /3 label the components of 
the spin operator. The subscript "0" means calculation 
in the case of vanishing interchain coupling and in the 
last relation the known result for the spin correlation, 
function of the homogenous Heisenberg chain is used.t3 
From conformal field theory it is known that this corre- 
lation function decays algebraically like 1/r^^, where x is 
the scaling dimension when we disregard multiplicative 
logarithmic corrections. We therefore conclude that the 
scaling dimension of the interchain coupling is x = 1 and 
represents a relevant perturbation of the critical system. 
From scaling relations we find, again ignoring logarithmic 
corrections, the ground-state energy of this system 
behaving like 



Eo cx J^- 

and a gap A is opening with 



Acx |Jj,| — 



Jt. 



(5) 



\Jy\ , (6) 



where d = 1 + 1 is the dimension of the correspond- 
ing classical model. The existence of a gap_£or the two- 
leg ladder has also been shown numerically.E^ In general. 



N/2 

h = jY,sI,{sI,_, 



^i) 



(7) 



r=l 



and by calculating again the energy-energy correlation 
function, we find that this perturbation also has scaling 
dimension a; = 1. That leads to the same conclusions as 
in the Ja;/Jj,-model. 

The situation is much more complicated in the meander- 
model, because there is not only an interchain cou- 
pling, but also a coupling between third-nearest neigh- 
bors within every chain as shown in figure ^. First, we 



FIG. 6: Two smooth chains from the meander-model 



want to investigate the intrachain coupling. The operator 
of this perturbation is given by 



N/2 

H = S2rS2r+3 



(8) 



and we can calculate the corresponding energy-energy 
correlation function 



^ 2r ' 



(9) 



where the value x = 1/2 of the scaling dimension of 



the singlet operator S2rS2r+3 has been employed.^ This 



nsipn 01 
LBThis 
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means that the NNNN couphng is relevant. By simply 
applying the scaling relations, we conclude that it de- 
stroys criticality and a gap opens with A cx | Jp^'^. But 
there might be some doubt if this scenario is correct, 
because if we suggest a short range Neel order on the 
critical chain the NNNN coupling is not frustrating. On 
the other hand, if the NNNN coupling is as strong as 
the NN coupling, this chain is equivalent to a 2-leg lad- 
der, which does show a gap. We therefore have used the 
transfer-matrix DMRG and the standard DMRG to test 
numerically the predictions from scaling. Before we enter 
the numerical part, we have to analyse the other pertur- 
bation in the meander-model caused by the interchain 
coupling. It turns out that this is again relevant with a 
scaling dimension x = 1 as in the other two models. 
The conclusion from scaling arguments is therefore that 
the 3x1 Jy- and the stair-model show one-dimensional be- 
haviour only at the decoupling point. Because in the 
meander-model the intrachain is more relevant than the 
interchain coupling, the scaling arguments suggest the 
existence of a disordered phase between the decoupling 
point and the phase with two-dimensional antiferromag- 
netic long-range order. This will be further investigated 



in sections HI and IV 



III. NUMERICAL INVESTIGATIONS 

To prove the scaling argument for the NNNN coupling 
in the meander-model, we have used two numerical meth- 
ods. The first one is the so called transfer- matrix pMRG 
(TMRG), which combines White's DMRG ideakl with 
the transfer-matrix approach.tS This methcjdj|has .been 
applied to different quantum chains beforellZniSllj and 
yields very accurate results for finite temperature. It 
is particulary suited, because the thermodynamic limit 
in quantum space can be performed exactly. Before we 
state the results, we write down the considered Hamilto- 
nian explicitly: 



N 



N/2 



2r+3 



(10) 



The relevant case in this context is a NNNN coupling J2, 
which is much smaller than the NN coupling Ji. 
The free energy as calculated with the TMRG for 
J2I J\ = 0.25 is shown in figure 0. We can now determine 
if there is a gap or not, because we know from scaling re- 
lations that in the low-temperature limit the free energy 
of a gapless, critical system scales like f{T) cx eo — a • 
with eo being the ground-state energy, whereas f{T) cc 
eo — a ■ r3/2g-A/T j£ there is a gap. We tried to fit a 
quadratic function to the data and noted that this is im- 
possible, whereas a function as expected for the gapped 
case fits perfectly with values eo = -0.46873 ± 0.00002, 
a = 0.29 ± 0.03, A = 0.23 ± 0.02 and errors, which are 
determined by a variation of the fit-region (see inset of 

fig- §■ 




FIG. 7: Free energy for the Heisenberg chain with NNNN 
coupling calculated by TMRG with J2/J1 = 0.25, m = 40 
states kept in the DMRG and a Trotter parameter of e = 
l3/M = 0.05. The inset shows a low-temperature fit with 
/ oc eo - aT^/^e"'^/^ and T G [0, 0.2]. 



This means that at this strength the NNNN coupling has 
really destroyed criticality. To test the scaling argument 
further, we also applied a standard DMRG program to 
this problem. When using the same parameters J2/J1 = 
0.25, we find a gap depending on the length of the chain 
as shown in figure ^ An extrapolation L — > 00 then leads 
to a gap ApBc ~ 0.23652 ± 0.00064 and a ground-state 
energy epBc = —0.46841 ± 0.00816 if periodic boundary 
conditions (PBC) are applied. For open boundary con- 
ditions (OBC), we find a gap Aqbc = 0.23834 ± 0.00027 
and a ground-state energy eoBC = —0.46843 ± 0.00003. 
Consequently, there is a good agreement between the 
numerical results from the two different methods. Up 
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FIG. 8: DMRG calculation of the gap AE for finite chains 
with length L for open and periodic boundary conditions us- 
ing an extrapolation in the number of states m. The lines are 
guides to the eye. 

to now, we have only stated that the system has a gap 



for one special choice of parameters. Because the infi- 
nite chain with J2 = is a critical system, we can use 
the RG to study the behaviour of the free energy de- 
pending on the length of the chain and of the NNNN 
coupling J2. In general, we can linearize the RG trans- 
formation in the vicinity of a critical Hamiltonian, which 
is a fixed point of the RG flow, and find that the free 
energy per lattice site for a classical system scales like 
/(51, ■■■:9n) = b-'^fib'^^gi, ■ ■ ■ , 6'^"g„) if the RG is ap- 
plied 1-times. Here gi denotes a linear scaling field, is 
the eigenvalue of the RG transformation and b is the scal- 
ing factor. Because the quantum chain with length L is 
equivalent to a classical system with volume V = L x (3, 
with P being the inverse temperature, the dimension d 
is equal to 2 and the relevant scaling fields at T = are 
1/L and J2. It follows that 

fiL,j,) = b-'^f(^b^^'^y^'j2^ (11) 

and by choosing b^^'^ = L we get 

At J2 = this reduces to f{L,0) = const • L~'^/^^ and 
because the ground-state energy per lattice site scales 
like we conclude that Ai = 1. 
By setting 6^^' = Jj"^ we find the relation 

/ (L, J2) = jf 'V (l- V2-'^/^^ 1) . (13) 

When inserting the known result A2 = 3/2, we can state 
that there must exist an universal scaling function $ with 

J-'/'-fiL,J,) = ^[L-^J-'^') . (14) 

We can do similar calculations for the gap A and get from 
A (L, J2) = 6"'A the scaling relation 

J2"'/^.A(L,J2) = |.(l. J2'/^). (15) 

To test this, we have applied the DMRG to 862 different 
chains with lengths up to 122 sites and J2 G [0.00005, 0.8] 
(see figure ^). By these calculations the scaling relation 
( p^ ) is confirmed in a convincing way. Note that practi- 
cally all data points collapse on a one-dimensional man- 
ifold. Some minor deviations are noticeable and can be 
explained by higher order terms in the finite size scaling 
and an effective exponent 0.662 instead of 2/3 (see in- 
set of fig. H) . It is also possible to determine two scaling 
regions. For large lengths L of the chain the plotted func- 
tion saturates, indicating that the relation AE cx J2^^, 

derived from scaling arguments for the infinite chain, re- 

2/3 

ally holds. For small L • there is a linear regime, 
showing that the finite size gap proportional to 1/L is 
the dominant contribution in this region. 
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FIG. 9: DMRG results for 862 chains with different length L 
and NNNN coupling J2. Two scaling regimes are visible as 
discussed in the main body. The inset shows the gap extrap- 
olated to the thermodynamic limit versus coupling J2. 

IV. MAPPING ONTO A NL(J-MODEL 

In this section, we want to discuss the possible tran- 
sition from the magnetically ordered phase to a disor- 
dered phase driven by the dimerization S by using the 
nlcr-model as a low energy effective theory. The easiest 
way to get a path integraL-f pr ,, the . considered models is 
the use of coherent states ri]'EllE§E3 Spin coherent states 
\n) form an overcomplete basis set and are generated by 
a S'[/(2)-rotation of the highest weight state \s, s) 

\n) :=e''^("°^")'^|s,s) , (16) 

where no is an unit vector along the quantization axis 
and cos(0) — n ■ Uq. By using the Trotter formula and 
inserting the identity operator, the partition function can 
be written as Z = J Z3ne~'^^["I with an Euclidian action 
given by 

/■/? 

SE[n] = -isSwz + dt{n\H\n) . (17) 
Jo 

Swz is a topological term (Wess-Zumino term), which 
arises from Berry phases and can be expressed as 

Swz[n] = 'y^ dr dt n ■ (dtn x dru) (18) 
r Jo Jo 

with the boundary conditions n(t, 0) = n(i), n{t,l) = 
no and n(0,T) = n{(3,T). 

A. Jx / Jy-moAe\ 

With flQ being the lattice constant, a site on the lat- 
tice G can be described by i = o + aoX]a=i*a^" with 
spatial unit vectors e°. Using this notation, which is 
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more suitable for the following calculations than , the 
Hamiltonian of the model can be represented as 



H = ^{ J^S'(i)5'(i + aoe^) + JyS{i)S{i + aoe^)} 



(19) 



We assume periodic boundary conditions and an even 
number of sites N in each direction. By defining R — 
Jy/Jx and using the coherent state relation (n| S \n) — 
sn, the second term in (|l^) can now be easily calculated 
leading to the Euclidian action 



Ssin] = ~isSwz[n 
r/3 



(20) 



f 

+Jxs'^^^ / dt {n{i)n{i + aoe'') + Rn{i)n{i + aoe^)} 



Now we use the well-known ansat: 



.atJi 



n{i) ^ p{i)m{i)^l - al'^l^ii) + (21) 



with 



p{i) = {-1)^1^1^' 



(22) 



taking into account the short-range Neel order due to 
the antiferromagnetic exchange. Here m is the order pa- 
rameter field and I represents the rapidly varying but 
small part. The constraint = 1 leads to = 1 and 
m ■ I = and we will expand ( pl| ) up to quadratic order. 
Because the Wess-Zumino term is independent of the mi- 
croscopic details of the spin model, we want to discuss 
this term more generally in d dimensions. Starting with 
( |l8| ) and using the ansatz (^ leads to 

isSivzin] ~ is / d'^x / dtl-{mxdtm) (23) 
Jv Jo 



N 



i2,--- ,id = i 



where 9 — 27rs and k{i2, ■ ■ ■ ,id) is the winding number 
or Pontryagin index of the field m defined by 

ki^2, • • • = ^ J dxi j dt dx^m ■ {m x dtm). (24) 



In one spatial dimension the second term in (^3|) is re- 
sponsible for the different physics of chains with inte- 
ger and half-integer spin. If the m-field is smooth, the 
integer-valued Pontryagin index k(i2,--- ,id) must be 
a constant rajid hence this term cancels out in higher 
dimensions.EJ Note that we have to treat this term more 
carefully for the anisotropic models considered here, be- 
cause the m-field may no longer be smooth in each di- 
rection. Using the same ansatz ( pl| ) for the part of the 
action depending on the Hamiltonian of the system and 
integrating out the i-field, results in an effective action 
for the low lying excitations 



S\m] 



[dxmf + R{dymf + \{dtm) 



(25) 



with a spin stiffness ps = JxS^, a transversal magnetic 
susceptibility xi. = '^Jx<io{^-\- R) and a spinwave velocity 
Vs = y/psXT- We now rescale the imaginary time by 
xq = Vs ■ t leading to 



Sniaim] 



1 



d''x{idxomf+ {dxmf+ Ridymf} , 

(26) 



where the dimensionless coupling gQ is defined by 

2 , 

50 = -Vi + R ■ 

s 



(27) 



We want to discuss ( |2^ ) fallowing some arguments given 
by Affleck and Halperin.tJ Because the coupling in y- 
direction may be arbitrarily weak, a continuum repre- 
sentation may not be justified and we therefore rewrite 
the action as 



R 

d^x <j (dxaTrinf + {dxTUnf + — (m„+i 



(28) 



After a rescaling y' ~ y/\/R the momentum space UV- the other directions. In a Wilsonian RG step, where 
cutoff in y-direction is now smaller than the cutoff in the higher momentum modes are integrated out, only 
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ko and contribute, so that the RG is essentially two- 
dimensional. Only if the momentum scale has been low- 
ered so that also ky has components in the shell, we have 
to switch to the three-dimensional RG. During the two- 
dimensional RG, we also have to consider the rescaling of 
the m-field which is given by m„ (A'/A)^m„. Here 
A is the UV-cutoff before renormalization and A' the re- 
duced one after a RG step, x is the scaling dimension of 
the m-field which is equal 1/2 for the Heisenberg model. 
The scale A', where we have to switch from two- to three- 
dimensional RG, is therefore given by the condition 



i?^(A7A)2-AV2go « A'V2g{A') 



(29) 



Here g(A') is the renormalized coupling constant when 
the momentum modes have been integrated out down to 
the cutoff A'. Note that the assumption that the in- 
field is smooth on the scale of the lattice spacing is no 
longer justified in y-direction. Instead of canceling out, 
the second part of (^3|) leads to an independent wind- 
ing number for each chain with topological angle 6 — n. 
Therefore g flows in two-dimensions to the marginally 
stable fixed point 52(0) of the s = 1/2 chain. As a con- 
sequence go/g{A') is always of order one and ( p9| ) sim- 
plifies to A' = RA. The coupling constant g{A') then 
acts as bare coupling constant for the three-dimensional 
RG fiow. Therefore the ground state is ordered if g{A') 
is smaller than the critical fixed point gc of the three- 
dimensional RG, whereas it is disordered if g(A') > gc- 
If even 172(0) < gc, the system is always ordered for a 
nonvanis hipg . R . Numerical calculations using different 
methodsllSE3c3 give strong evidence that this model or- 
ders for arbitrarily weak R, meaning in the language of 
RG that (72(0) seems to be smaller than gc (see the cor- 
responding fiow diagram fig.nffl). 



g(RA) 

— • — 



g,(0) 



Jx / Jy-'^nodel with the sum in (|22) now running up to d, 
the Hamiltonian can be expressed by 



= ^ [1 + pii)6] S{i)Sii + aoe") . 



(30) 



ieG a=l 



Using again relation (^^ and also the ansatz (|2l|), an 
Euclidian action for this model depending on the unit 
vector field m and the orthogonal vector field I is derived. 
By the same arguments given in the chapter before, the 
Wess-Zumino term vanishes in spatial dimensions greater 
than one under the assumption that m is smooth on 
the scale of the lattice spacing, but is important in one 
dimension. Again, we integrate out the rapidly varying 
but small Z-fields. The result is an effective action 



S[m] 



2 



r0 

(fx / dt 

V Jo 



(31) 



with a spin stiffness ps — Js (1 — ^), an induced 
anisotropy 7"'' = 5 



1 — 5 ), a transversal mag- 
netic susceptibility xa- — 4 Jdog and a spinwave velocity 
'Vs — ^/PsX±- But there arises also an imaginary con- 
tribution proportional to 6 from the Z-field integration, 
which can be expressed as 

—j- J dxi J dt dx^rn ■ [m x dtm) + ■ ■ ■ 
+ X] y"^^d j dt dx^m ■ {m X dtm) (32) 
= -J- ki\T'2, ■■■ H 



kdiii, ■ ■ ■ ,id-i) ■ 



1-1, ■■■ ,4d-l 



If the m-field is smooth, every ki must be a constant. 
Under this assumption the topological term simplifies to 



FIG. 10: RG flow for the / Jy-model 



B. stair- model 

We generalize the stair-model to d dimensions assum- 
ing a hypercubic, bipartite lattice G = A (B B with pe- 
riodic boundary conditions and an even number N of 
sites in each spatial direction. We do this, because the 
model reduces in one dimension to a dimerized chain 
and we want to compare the results, especially the topo- 
logical terms, for the dimerized chain and for the two- 
dimensional model. Using the same notation as for the 



top 



—N" 
d 



(33) 



In one dimension, where the stair-model corresponds to 
the dimerized chain, we find a total topological contribfc 
tion of Stop = iS{l + 6)k as has been calculated before.EZi 
Before we discuss the two topological terms (|2^ ) and ( |3^ ) 
in two dimensions, we look at the additional anisotropy 
in the action (31) expressed by the matrix 7°''. This ma- 
trix is symmetric and becomes the identity if 5 goes to 
zero. 

In two dimensions a 45°-rotation diagonalizes this matrix 
and with xq = Vg ■ t the action of the two dimensional 
model is given by 



d^x —{damf 



(34) 
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There are now different spin stiffnesses in the spatial di- 
rections given by 



p\ = Js^ = Ps 



2-6^ 



; Ps= Ps 



2-2(5^ 
2 -,52 



(35) 



and = Ps- By a rescaling of the imaginary time. ( |3^ ) 
is transformed into 

S[m] = fd^S:{(d,mf + (1 - 6^){dym)^ + {d,,mf} 

(36) 



with a bare coupHng given by 



So 



2V2 



(37) 



Again the coupHng in y direction may be arbitrarily 
weak, and we therefore have to use a discrete version 
of (^6|) like for the / Jy-mode\. As a consequence a 
two-dimensional RG has to be used until the UV-cutoff 
is lowered to A' = (1 — S'^)A. The coupling g{A') is then 
the bare coupling for the three-dimensional RG. Now, 
we have to remember that there are also two topological 
terms (^ and ^ present. Because we have stated that 
it is necessary to use a discrete representation instead of 
derivatives for the weak couplings, neither the winding 
number in x- nor the winding number in y-direction is 
well defined any longer. We therefore have performed 
an alternative mapping to a nlcr-model starting with a 
slightly modified version of this model, where the chains 
formed by strong bonds are smooth and along the x-axis 
and the weak bonds form zigzag-chains (see figure Iq). 
The result is again an anisotropic nlcr-model like ([26[), 
but now without a topological contribution proportional 
to S and with a winding number in the topological part 
of (HI) calculated along the strong bonds. The situa- 
tion is therefore exactly the same as in the J^;/ Jj,-model 
and if we accept 52(0) < gc as an universal property, we 
conclude that the ground state of this model is always 
antiferromagnetically ordered for S G [0,1), and this or- 
der only vanishes at 5 = 1, where the model consists of 
uncoupled critical chains. 



C. plaquette-model 

The model is described by the following Hamiltonian 

+ [l+Py{i)S]S,,jS,,j+i} , (38) 



where we have defined Px (i) = (—1)'=" andpj,(i) = (— 1)*«. 
Substituting (|38| ) into (j^), we derive again a path inte- 
gral formulation. By using the same ansatz as before, 
we notice that the dimerization S does only contribute 
in (daTTi)^ ■ I and higher orders. This is not surprising, 
because the relevant low energy modes are not only near 
wave vectors (tt, tt) and (0, 0) but also near (0, tt) and 
(tt, O). We therefore generalize ([21]) in the following way: 



n{i)=p{i)m{i)\Jl - al'^[lQ{i) + Px{i)lx{i) + Vy{i)ly{i)\^ 
+4[lo{i) + px{^)l (*)] (39) 

The constraint = 1 implies — 1 and m ■ li = 0, 
where li denotes any of the three fluctuation fields. By 
an expansion up to quadratic order in li, dam and dtm, 
we get again the effective action ( p5| ) but now R = \ and 
the other parameters are defined by ps = Js^{l — S"^), 
X± = 8 Jaq and Vs = y^pTxT- By rescaling the imaginary 
time axis this is transformed to the standard nla-model 

S[m] = / d^x{^^m){^^'m) (40) 

with a bare coupling 



90 



2V2 



1 



(41) 



All topological terms vanish for this model and - differ- 
ent from the two models considered before - the action 
is isotropic in the spatial directions. So the 3d RG has 
to be used from the beginning and whether or not the 
model has an ordered ground state depends on whether 
or not go < gc- Because go 00 for (5 — > 1, the or- 
der vanishes already before the plaquettes decouple. The 
precise value of Sc € (0, 1) cannot be determined within 
the RG treatment. 



D. meander- model 

The Hamiltonian looks quite similar to (|38|), but in- 
stead of an alternation with Px(i) in x-direction there is 
one with p(i). As a result the Zy-field is unnecessary and 
can be left out. By using this reduced form of (^9|) and 
integrating out the Iq- and Z^j-fields an effective theory is 
derived 
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i=l 



(42) 



(fx 
V Jo 



(l-^)(a.m)2 



{l-5^){dymf + \{dtm) 



where ps — Js^, x± — 8 Joq and Vg — •JpTxZ- Like in the 
stair-model there are different spin stiffnesses in the spa- 
tial directions and also a topological contribution propor- 
tional to (5 like (|3^), but in this model only the winding 
number in x-direction is involved. As in the stair-model, 
this winding number is not well defined. The same cal- 
culations within the modified model, where the chains 
formed by strong bonds are again smooth and along the 
X-axis (see figure show that the topological term pro- 
portional to (5 vanishes and the winding number in ( |2^ ) 
is calculated along the strong bonds. We therefore ig- 
nore the topological term in ( ^2| ) from now on and use 
a 2d RG flow with 9 = n later on. A rescaling of the 
imaginary time axis then leads to 



S\m] = - 



1 



2ao5o, 



2(52 



2-52 



{dymf 



where the bare coupling is given by 

2V2 I 2 



9o{S) = 



2 -(52 



-(atm)2l 
(43) 



(44) 



Because of the anisotropy in y-direction the RG is, like 
in the Jx/Jy- and stair-model, two-dimensional at the 
beginning. But here the bare coupling gQ{6) is increased 
with increasing 5. Therefore two scenarios are possible: 
If a (5c exists so that .go(<5) > .92(0) if (5 > 5c, the coupling 
is driven by the 2d RG flow towards inflnity and (?(A') 
with A' = (2-2(52)/(2-52)A, which is the bare coupling 
for the 3d RG, is then greater than gc- Therefore the 
ground state is disordered for 5 > 5c- If however even for 
5 ^ \, gQ{5) remains smaller than 52(0), the ground state 
is always ordered for a nonvanishing interchain coupling. 
In summary, we cannot decide within the RG treatment 
if an extended disordered phase exists or if there is long- 
range order for all 5 £ [0, 1) as in the Jx/Jy- and the 
stair-model. 



The second reason is that we are interested in the ques- 
tion if these models can be the result of a dynamical pro- 
cess i.e. a spin-Peierls transition. Because the calculated 
ground-state energy for the isotropic two-dimensional 
Heisenberg antiferromagnet of Eq/NJ k, — 0.65 7|9 devi- 
ates less than 3% from the best numerical resultsOa and 
also the sublattice magnetization m — 0.3034 agrees very 
well, we expect that LSWT gives reliable results near 
the isotropic point. On the other hand, the results for 
large dimerizations have to be regarded with care, be- 
cause LSWT fails in one dimension. |— ■ 
We use the Holstein-Primakoff transformationEa to map 
the spin operators onto Bose operators: 



a a, 



2s 



1 



2s 



(45) 



By expanding the square root in 1 /s and taking only the 
lowest order into account, LSWT is reached. 



A. stair-model 

As before, we want to consider the stair-model gener- 
alized to d dimensions. Because the lattice G = A(B B is 
bipartite, it is possible to write the Hamiltonian as 

d 

H = ^^[J(l + <5)S'(7;)S'(i + aoe") 

+ J{l-5)S{i)Sii-aQe^)] , (46) 

which is useful for LSWT. Starting point is again a 
Neel ordered state, and we therefore apply independent 
Holstein-Primakoff transformations to the two sublat- 
tices A, B: 

X e A: S'-{x) = s - a+{x)a{x) ; {x) = V2sa+{x) 



(47a) 



V. LINEAR SPINWAVE THEORY (LSWT) 

In the section before, we have mentioned that it is im- 
possible to determine a value for the critical coupling gc, 
or equivalent a value for the criticial dimerization 5c, for 
the plaquette- and the meander- model from RG. This is 
one reason to apply LSWT onto the considered models. 



xe B : S''{x) = -s + b+{x)b{x) ; [x) = V2^b{x) 

(47b) 

Taking the Fourier transform the Hamiltonian is bilinear 
and given by 

H ^ -NJs^d (48) 
-I- 2Jsd"^{a+ak + b+bk + Aakbk + A*a+b+} , 
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where 

d 2^r. 



the definitions = 
_j^sin(fcmao) ^J^d A — 



- aEf=iCos(fc,ao), (3k = 
"fk + ^Pk have been used. 



By means of a Bogoliubov transformation this is easily 
diagonahzed leading to 



H 



-NJds{s 
2Jsd'Y^ 



1) 



(49) 



^-ill-5'Pl){clck + dldk + l) 



1 Sublattice magnetization 

The sublattice magnetization m is given by 
N 



m 



^{alak) 

k 

_ N 1 \p / 1 



(50) 



By replacing the sum by an integral this can be evaluated 
in principle in any dimension d. However in one dimen- 
sion this integral is divergent. This is not astonishing 
due to Coleman's theoremE^ stating that the continuous 
S'C/(2)-symmetry cannot be broken in one dimension. 
In two dimensions the dimerization 5 reduces the magne- 
tization as expected and aX 5 — \ there is again the one 
dimensional infrared divergence. Nevertheless, we can 
take the value of the dimerization where m vanishes as 
indication for the breakdown of Neel order. For the stair- 
model we ggt = 0.8286 or i?c = (1 - 5) /{I + &) = 0.094 
(see figure |lj) . 

To check the validity of the LSWT, we want to use an 
argument given by Sakai and Takahashi.E^ In spinwave 
theory the number Ux of bosons per lattice site is not 
restricted. That means that there are unphysical states 
in this theory, because in the original spin system the 
condition Ux < 2S holds. A possible estimation for the 
validity of LSWT may therefore be given by 



+ Aux <2s ^ (Sx) - Ahx > 



(51) 



This predicts that LSWT is valid for this model from 
(5 = up to (S « 0.65 (see figure |ll|). Therefore the 
calculated Sc is in a region, where LSWT is not reliable. 
We can only conclude that Sc must be greater than 0.65. 



2 Distortion due to spin-phonon coupling 

The dimerization leads to a gain in magnetic energy on 
the one hand, on the other hand it costs elastic energy. To 
answer the question if a spin-Peierls transition towards 
a distorted lattice is possible, we expand the calculated 
ground-state energy up to quadratic order. Because the 
linear term vanishes for symmetry reasons, we get 



0.4 



0.2 - 



0.1 



-0.2 



-0.4 




is:) 

\ + {SD - An, 



FIG. 11: Sublattice magnetization and condition (^l|) for the 
validity of the LSWT, both for the stair-model. 



with eo ~ —0.6579 and A = 0.261. As mentioned in the 
introduction (|^) the displacement of an atom is propor- 
tional to 5 and the mean-field treatment of the phonons 
leads to a classical elastic energy proportional to 5^. 
Therefore the lattice is only distorted if 



m J^Wj^^^^^TT, tt) 
2" A2_(^,7r) 



< 0.261 , 



(53) 



whereas the lattice is unchanged otherwise. Here m is the 
mass of the moved atom and lo the phonon frequency of 
the responsible phonon mode. If this condition is fulfilled, 
there is a coexistence of spin-Peierls dimerization and 
Neel-like long-range order. The situation here is quite 
different from that in one dimension. For the dimerized 
chain the magnetic energy scales like Emag <5'*^^ and 
because the clastic energy is proportional to 5^, there is 
always a distortion. In the two-dimensional stair-model 
the magnetic and elastic energy behave like 5'^ and there- 
fore a distortion of the lattice is only possible if certain 
conditions are fulfilled. Because we expand around the 
isotropic point where LSWT results agree well with nu- 
merical results, the quantitative results for the condition 
of dimerization are expected to be reliable. 



B. Ja;/ Jy-model 

Because this model is j«ell known and LSWT has al- 
ready been applied to it,L3 we only want to briefly state 
the results. The bosonic Hamiltonian is again easily di- 
agonahzed by a Bogoliubov transformation and can be 
written as 

H = -N{Jx + Jy)s{s + 1) (54) 
+ 2sJx{l + R)'^-fk{ctck + d+dk + l) 



with 7fc 



{cos{kx 



- Rcoa{ky)y 



and R - 



(52) 
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1 Sublattice magnetization 

If we again calculate the sublattice magnetization and 
take its vanishing as indication for the breakdown of Neel 
order, we find a critical Rc — 0.0337. We also proved 
by using condition ( pT] ) the reliability of the spinwave 
approach. It seems to be valid if i? > 0.1 and therefore we 
only can conclude that the critical coupling Rc must be 
smaller than 0.1. As mentioned before, numerical studies 
such as series expansionalJ and Monte Carlo studiesE£l 
seem to indicate that Rc is equal to zero. 



2 Distortion due to spin-phonon coupling 

By using the relation S = (1 — i?)/(l + i?) we switch 
again to the description of the model with the parameter 
S. In the expansion of the ground-state energy (|2|) the 
parameter A is now given hy A = 0.146. A spin-Peierls 
transition is only possible if the gain of magnetic energy 
is larger than the cost of elastic energy leading to the 
condition 



(2cL,(l, 0) + cl_{l, 1) - cljl, 1)) < 0.146 



(55) 



where c(l, 0) and c(l, 1) are the phonon velocities in di- 
rection (1,0) and (1,1). 



the Fourier transform the Hamiltonian with additional 
field is given by 

H = -2NJs^ - NB^s 

+ 4Js^ {ofcOfc + b+bk + c^Ck + d^dk 

k 

+ Al [ofcCfc + b+d+] +h.c.+Al [akdk + b+c+] +h.c. 
+ §f; H^k + b+bk + c+ck + d+dk] } , (56) 

where 7^ = i cos(A;aoe°), /3£ — ^sin(A;aoe") and A^. = 
lk + 5(3l. 

This has to be diagonalized under the subcondition that 
the new operatacs fulfill again Bose commutation rela- 
tions leading toEll 



and 



(I) T-^UJT ^UdJ 



(II) T+JTJ = 1 



(57a) 



(57b) 



where T-Lu denotes the diagonalized Hamilton matrix and 
J is given by 



J = 



-1 
0-100 
10 
1 



(57c) 



C. plaquette- model 

The situation is more complicated for this model, 
because the unit cell includes four sites. Therefore we 
have to introduce four kinds of bosons. This happens as 
follows 

A: r e (2i, 2j); S"^ = s - a+a^; S- = ^/2's a+ 

B: r e (2i + l,2j + l); 5" = s - ^ ^/2'sb+ 

C: r e {2i, 2j + 1); 5^ = ~s + c+Cr] = y/Ys Cr 

D: r e {2i -f 1, 2j); 5^ = -s + d+dr\ S- = ^2^ dr , 

where A, B, C, D enumerate the four sublattices. 
In principle it is possible to diagonalize every Hamilto- 
nian of an assembly of N bilinearly intera|Cting bosons 
or fermions what is well-known since long.cJ But espe- 
cially for bosons it is often complicated to construct the 
transformation matrix T between new and old operators, 
because the transformation is not unitary. We there- 
fore add an "antiferromagnetic field" B\, which allows 
us to calculate the sublattice magnetization without do- 
ing this canonical transformation explicitly. After taking 



From equation (II) it follows that T is an element of 
the pseudounitary group ?7(2, 2), whereas (I) implies the 
secular equation 



det{nJ ~Xil) = 



(58) 



for the diagonalization problem. The eigenvalues of Ti. are 
then given by |A,i|. By solving the secular equation ( [S^ ) 
we get the diagonalized Hamiltonian of the plaquette- 
model in LSWT: 

H = -2NJs{s + l) - NB%{s + ^) 
+ 4Js^{Ai [A+Ak + C+Ck + 1] 

k 

+ A| [B+Bk+V+Vk + 1] } (59a) 



= V(l + Si)2-(|A-| + |A^,|)2 



X2 



a + B%r-i\Al\-\Al\r (59b) 



with the new Bose operators Ak, Bk, Ck, T)k and Bj^ ~ 
B\lAJs. 



1 Sublattice magnetization 



The ground-state sublattice ma gnet ization per lattice 
site is now easily calculated from (|59a|) by the derivative 
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Ib^ -0' Setting s = 1/2 and replacing the sum 
through an integral the sublattice magnetization is given 

by 



Tr/2 /.■rr/2 



-7r/2 J-w/2 



(60) 



where AJ, — \\{B\ — 0). This decreases with increasing 
5 and vanishes at 5c — 0.798 or aquivalently Rc — 0.112. 
This value agrees with-ihat given in a paper by Koga, Ku- 
mada and KawakamiE3 The same authors have also used 
a series expansions starting from uncoupled plaquettes 
to determine the critical coupling and get 5c ~ 0.3. This 
value is in good agreement with results from Monte-Carlo 
calculations, showing again that the results in LSWT at 
high dimerizations have to be considered with care, be- 
cause of the unphysical states in this approach. 



2 Distortion due to spin-phonon coupling 



We set B\ = and expand again the ground-state en- 
ergy up to quadratic order. For this model the parameter 
A is equal 0.174. Therefore the lattice is distorted if the 
condition 



AL.(^,0) ^ A2 (0,7r) 



< 0.174 



(61) 



is fulfilled. 



D. meander-model 

The sublattice structure of this model is similar to that 
of the plaquette-model and it is again necessary to intro- 
duce four kinds of bosons. With the same definitions as 
before the Hamiltonian of the model is given by 

H = -2NJs^ - NB\s 

+ 4,Js^^{a^ak + b^bk + c^Ck + d^dk 

k 

+ Al [flfcCfc + bkdk] + h.c. + Al [akdk + b'^c^] + h.c. 



+ ^ [a^ak + b+bk + c+Ck + d+dk] }■ 



(62) 



By solving the secular equation (|58|) this is diagonalized 
leading to a Ha milto nian in the new Bose operators Ak, 
Bk, Ck, T^k like (59a) but now with eigenvalues 



= \/ 1 + Si - \Al\' \Al\' 2 \Al\ lk{Al) (63) 



>^i=\ [^+B\) - \Ai\' - m + 2 \Ai\ mi) 



1 Sublattice magnetization 

Analogous to the plaquette model the sublattice mag- 
netization can be calculated by a derivative. The value 
for the critical dimerization determined by this calcula- 
tion is 5c = 0.898 or Rc = 0.054. This is like in all the 
other models in a region where the spinwave approach is 
no longer justified. 



2 Distortion due to spin-phonon coupling 

The parameter in the expansion is given hy A = 0.160 
and therefore the condition for the spin-Peierls transition 



m / J^a;^(7r, tt) 
T A2(7r,7r) 



AL.(0,7r) 



< 0.160 , 



(64) 



where we have assumed that the phonon frequencies and 
spin-phonon coupling constants for the longitudinal and 
transversal (tt, 7r)-phonon are identical. 



VI. CONCLUSIONS 

By using simple scaling arguments we have shown that 
for the Ja;/ Jj,-model and the stair-model, which consist of 
weakly coupled Heisenberg chains for large dimerizations, 
the interchain coupling is a relevant operator. Therefore 
only &i 5 = \ the ground state is of one-dimensional na- 
ture. Also in the meander-model the interchain coupling 
is relevant, but there is also a relevant intrachain cou- 
pling. By numerical calculations (DMRG + TMRG) we 
have shown in this more complicated case that the scal- 
ing behaviour for the intrachain coupling is predicted cor- 
rectly by the simple scaling argument. We also showed 
that it is possible to map all four models onto a nla- 
model. By applying the known results for the RG fiow 
in two and three dimensions we have concluded that the 
ground state of the Jx/ Jy- and the stair-model is Neel- 
like ordered for all 5 E [0, 1), whereas an extended dis- 
ordered phase exists for the plaquette-model. However, 
it was not possible to determine the critical dimerization 
5c G (0, 1) for the plaquette-model within the RG treat- 
ment. For the meander-model the RG gave no unique 
result. A phase diagram like for the Jx/ Jy- and stair- 
model, but also an extended disordered phase as in the 
plaquette-model are both possible scenarios. The second 
possibility seems to be more probable, as the intrachain 
coupling has scaling dimension a; = 1/2 and is therefore 
more relevant than the interchain coupling with a: = 1. 
If we start from 5=1 and reduce the dimerization, we 
might expect that there arises first a system consisting 
of weakly coupled 2-leg ladders. Because a 2-leg ladder 
with s = 1/2 is a gapped system, a small coupling be- 
tween these ladders can be treated within normal pertur- 
bation theory and does not change the global properties 
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drastically. At lower dimerization the gap closes and the 
system orders antiferromagnetically. If this picture is cor- 
rect, Sc must be smaller than one. 

To investigate this further and to determine a value for 
the critical dimerizations Sc, we also applied LSWT to 
the models. The following values for Sc are derived: 



model 


Sc 


Jx/ Jj^-model 


0.935 


meander- model 


0.898 


stair-model 


0.829 


plaquette-model 


0.798 



All these values have to be regarded with great care, be- 
cause LSWT allows also unphysical states and a simple 
argument has shown that there is a large contribution of 
these states at such high dimerizations. 
This problem does not occur near S = 0, where LSWT 
gives very precise results. We therefore believe that the 
conditions for the spin-Peierls transition, which we have 
obtained within LSWT, are qualitatively and quantita- 
tively useful. For all models we have found that elastic 
and magnetic energies scale like S^ at small dimerizations. 
A phase transition leading to a coexistence of spin-Peierls 
dimerization and antifcrromagnetic long-range order is 
therefore only possible for certain values of the micro- 
scopic coupling constants. From an expansion of the 
magnetic ground-state energy E{S) = eg — A ■ S'^ we have 
got the following values for the parameter A: 



model 


A 


stair-model 


0.261 


plaquette-model 


0.174 


meander-model 


0.160 


Jx / Jy-model 


0.146 



The gain of magnetic energy is therefore largest for a 
stair-like distortion of the lattice, what is caused by a 
transversal (tt, 7r)-phonon. What kind of distortion is en- 
ergetically prefered, depends also on the elastic energy, 
which in general is different for each model. However, if 
we assume that the elastic energy is equal for all models, 
we would conclude that the stair-like distortion is ener- 
getically prefered. Titis is in contradiction to a result 
by Tang and Hirsch,a who have studied the plaquette- 
and stair-like distortion by an exact diagonalization of a 
4 X 4-lattice and conclude by using the same assumption 
that the plaquette structure is prefered. However, the 
lattice they have considered is very small and they have 
not done any finite size scaling so that we believe our 
result is more reliable for the infinite lattice. 
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